Vector meson form factors and their quark-mass dependence 
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The electromagnetic form factors of vector mesons are calculated in an explicitly Poincare co- 
variant formulation, based on the Dyson-Schwinger equations of QCD, that respects electromag- 
netic current conservation, and unambiguously incorporates effects from vector meson poles in the 
quark-photon vertex. This method incorporates a 2-parameter effective interaction, where the pa- 
rameters are constrained by the experimental values of chiral condensate and f^. This approach 
has successfully described a large amount of light-quark meson experimental data, e.g. ground state 
pseudoscalar masses and their electromagnetic form factors; ground state vector meson masses and 
strong and electroweak decays. Here we apply it to predict the electromagnetic properties of vector 
mesons. The results for the static properties of the p-meson are: charge radius (r 2 p ) = 0.54 fm 2 , 
magnetic moment p, — 2.01, and quadrupole moment Q = —0.41. We investigate the quark mass 
dependence of these static properties and find that our results at the charm quark mass are in 
agreement with recent lattice simulations. The charge radius decreases with increasing quark mass, 
but the magnetic moment is almost independent of the quark mass. 

PACS numbers: ll.10.St, 13.40.Gp, 14.40.-n 
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I. INTRODUCTION 

Hadron form factors provide an important tool for un- 
derstanding the structure of bound states in QCD. The 
coupling of a (virtual) photon to a composite particle de- 
pends on its internal structure. Even static properties 
such as the charge radius and magnetic moment are sen- 
sitive to the underlying QCD dynamics. Thus it is not 
surprising that there have been numerous studies of the 
electromagnetic form factors of the nucleon, both the- 
oretically and experimentally. Also the form factors of 
pseudoscalar mesons, in particular of the pion, have been 
studied extensively [J 0,11, II S, !, 0, 0, St E3, [H • There 
are fewer studies of the vector meson form factors [l2l . [T3 , 
0, [H[ , though in recent years there has been a renewed 
interest in these form factors [H, Q3, [U, EI S3, S3 ■ How- 
ever, the results of these theoretical studies appear to 
suffer from a rather large model dependence, e.g. the re- 
sults for the quadrupole moment Gq (0) differ by a factor 
of two between different theoretical calculations, and the 
situation gets worse as one moves away from Q 2 = [T3| • 

Here, we calculate the electromagnetic form factors of 
the p and both the neutral and charged K* mesons us- 
ing the model proposed in Ref. [2^] . The parameters of 
this model were adjusted to reproduce the experimen- 
tal values for the chiral condensate, the pion mass and 
decay constant, and the kaon mass; calculations of other 
mesonic observables are predictions. The obtained vector 
meson masses and decay constants are in agreement with 
data [22[ as are the results for the strong decays of the 
vector mesons Sjj. The pion and kaon electromagnetic 
form factors [g, \A are also in excellent agreement with 
experiments [l(J, [Dj, as are the transition form factors 
such as the p — > [24[. This proclivity of theoreti- 
cal results to experimental data provides good reason to 
expect that the model will accurately describe the elec- 
tromagnetic form factors of light vector mesons. 



In Sec. [II] we review the Dyson-Schwinger equations 
used to calculate the quark propagators and the meson 
Bethe-Salpeter amplitudes. Next, we discuss the trun- 
cation and the model for the effective interaction, and 
the results for the meson masses, decay constants, and 
their quark mass dependence. We explicitly demonstrate 
frame independence: our results for physical observables 
are independent of the total meson momentum. This is 
important since in form factor calculations at least one 
of the mesons is moving. In Sec. IIV1 we briefly discuss 
the general form of vector meson form factors and the 
method of calculation; an essential element of our calcu- 
lation is the treatment of the quark-photon vertex. We 
present our numerical results for the form factors of the 
p and of both the neutral and charged K* mesons in 
Sec. El as well as the mass dependence of the static elec- 
tromagnetic properties of equal-mass vector mesons; we 
also compare our results with other calculations and with 
recent lattice simulations. Finally, some of the details of 
our calculation are given in the appendix. 



II. DYSON-SCHWINGER EQUATIONS OF 
QCD 



The Dyson-Schwinger equations [DSEs] are the equa- 
tions of motion of a quantum field theory. They form 
an infinite hierarchy of coupled integral equations for the 
Green functions (n-point functions) of the theory. Bound 
states (mesons, baryons) appear as poles in the Green 
functions. Thus, a study of the poles in n-point func- 
tions using the set of DSEs will tell us about properties 
of hadrons. For recent reviews on the DSEs and their use 
in hadron physics, see Refs. [H, S3, S3, SI- 
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A. Quark propagator 

The exact DSE for the quark propagator is 1 
Sip)- 1 = i jf Z 2 + m q (0 Z 4 + 

Zi [ g 2 D^(q) 7 ^S(k)ri{k,p), 



(1) 



where D fll ,(q — fc — p) is the renormalized dressed 
gluon propagator, and r*(fc,p) is the renormalized 
dressed quark-gluon vertex. The notation f, stands for 

J A c? 4 fc/(27r) 4 . For divergent integrals a translationally 
invariant regularization is necessary; the regularization 
scale A is to be removed at the end of all calculations, 
after renormalization, and will be suppressed henceforth. 
The solution of Eq. |T]) can be written as 



S(P) 



Z( P 2 



i p' + M(p 2 ) 



(2) 



renormalized according to S(p)~ 1 = if + m(£) at a suf- 
ficiently large spacelike £ 2 , with m(£) the current quark 
mass at the scale £. Both the propagator, S(p), and the 
vertex, r* , depend on the quark flavor, although we have 
not indicated this explicitly. The renormalization con- 
stants Z2 and Z4 depend on the renormalization point 
and on the regularization mass-scale, but not on flavor: 
in our analysis we employ a flavor-independent renormal- 
ization scheme. 



at discrete values P 2 = — M 2 of the total meson 4- 
momentum P. Here T is the Bethe-Salpeter amplitude 
[BSA] 

r(&o Ut , hn,P) = S(k out )' 1 x(kout,k i n;P)S(k in )~ 1 , (5) 

and the kernel K is the qq irreducible quark-antiquark 
scattering kernel. 

The meson BSA is normalized according to 

d 



2 P u = N c „ 
P dQ, 



k,q 



)x(<?out,<7in)] 



Tr[r(fci n , fc ut) S'(fcout) r(fc ou t, ki n 



Q=P 
(6) 



at P 2 = -M 2 , with fc out - k- m = Q = Qout - qin- The 
properly normalized BSA r(p ou t,Pin; P) [or equivalently, 
xiPa-atiPin', P)] completely describes the meson as a qq 
bound state. Mesons of different spins and parity are 
characterized by different Dirac structures, e.g. the BSA 
of massive vector mesons can be decomposed into eight 
Dirac structures u2fl 



r^(fc + ?7 P,fc-(l-77)P;P) = 

8 

^r.f(fc 2 ,fc-p ;?7 )T;(fc,p), 



(7) 



B. Mesons 

Bound states correspond to poles in ro-point functions: 
for example a meson appears as a pole in the 2-quark, 
2-antiquark Green function = (0|<7i<72<Zi?2|0). In 

the vicinity of a meson, i.e. in the neighborhood of 
P 2 = — M 2 with M being the meson mass, such a Green 
function behaves like 



G (4) 



X(jW,Pin;P) X(fcin, fcout; f) 
P 2 + M 2 



(3) 



where P is the total 4-momentum of the meson, p out 
and pi n are the 4-momenta of the outgoing quark and 
incoming quark respectively, and similarly for fc; n and 
fc out . Momentum conservation relates these momenta: 

Pout Pin — P — fcout kin ■ 

The function x(Pout,Pin! P) describes the coupling of 
the bound state to a dressed quark and antiquark. It 
satisfies the homogeneous Bcthc-Salpeter equation [BSE] 



r (pout, Pin, P) = 

/ K(p out ,pi n ; 

Jk 



"-out j "-in )x(fcout,fci n ;P), (4) 



1 We use Euclidean metric {Tu,')v} = 25 ^ 
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where T^(fc,P) are eight independent transverse Dirac 
tensors. The invariant amplitudes p are Lorentz scalar 
functions of fc 2 and fc • P, and depend on the momentum 
partitioning parameter rj. Physical observables however 
are independent of 77. 



C. Rainbow- ladder truncation 

A viable truncation of the infinite set of DSEs has 
to respect relevant (global) symmetries of QCD such as 
chiral symmetry, Poincare covariance, and renormaliza- 
tion group invariance. For electromagnetic interactions 
the truncation should also respect current conservation. 
These properties are built into the rainbow-ladder trun- 
cation [|| @, Ht| [13, HH . In this scheme, the kernel K 
of the meson BSE is replaced by an (effective) one-gluon 
exchange 



K{p ou t, Pin] fcout, fcin) 



(8) 



where q = p ou t — fcout = Pin — fcin, and a(q 2 ) is an effective 
running coupling; D { ™ e (q) is the free gluon propagator; 
and we choose to work in Landau gauge. The correspond- 
ing rainbow truncation of the quark DSE is 



Z ig 2 D^(q)Tl(k,p) 



(9) 
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This truncation is the first term in a systematic ex- 
pansion [29| of the quark-antiquark scattering kernel K; 
asymptotically, it reduces to leading-order perturbation 
theory. Furthermore, these two truncations are mutually 
consistent in the sense that the combination produces 
vector and axial-vector vertices satisfying their respec- 
tive Ward identities. In the axial case, this ensures that 
in the chiral limit the ground state pseudoscalar mesons 
are the massless Goldstone bosons associated with chiral 
symmetry breaking [3(1 [HI]. ^ n tnc vector case, this en- 
sures, in combination with impulse approximation, elec- 
tromagnetic current conservation 0, @] • 



III. MODEL CALCULATIONS 

The ultraviolet behavior of the effective running cou- 
pling is dictated by the one-loop renormalization group 
equation; the infrared behavior of the effective inter- 
action is modeled, and constrained by phenomenology. 
Here, we employ the model of Ref. [22| for a(q 2 ) 



A-Ka(q 2 



4lT 2 Dk 2 



-k 2 /u 



An 2 lm T{k 2 ) 



I In 



+ (i + k 2 /A 2 QCD y 



(10) 



and 



with F(s) = (1 - exp ^)/s, 7m = 12/(33 - 27V/), 

fixed parameters m t = 0.5 GeV, t = e 2 — 1, Nf = 4, and 
Aqcd = 0.234 GeV. The remaining parameters u> and D 
were fitted in Ref. L 22] to reproduce f n and the chiral 
condensate: ui = OA GeV and D = 0.93 GeV 2 . 



A. Results for light quarks 

With this model, we obtain good agreement with the 
experimental values for the light pseudoscalar and vec- 
tor meson masses and leptonic decay constants, see Ta- 
ble |TJ The current quark masses m u / d — 3.7 MeV and 
rn s = 83.8 MeV at the renormalization point C = 19 GeV 
were fitted [13] to the pion and kaon mass respectively. 
Using the one-loop expression to evolve these masses 
down to C = 2 GeV gives m u/d (l GeV) = 5.0 MeV and 
m s (2 GeV) = 118 MeV. 

These results show little sensitivity to variations in the 
model parameters (36| , as long as the integrated strength 
of the effective interaction is strong enough to generate 
an acceptable amount of chiral symmetry breaking, as 
indicated by the chiral condensate. This is not true for 
heavier states consisting of light quarks: e.g. the radially 
excited pion is quite sensitive to details of the interac- 
tion [13]. 

Not only the meson masses and leptonic decay con- 
stants, but also a wide range of other observables agree 
with experiments, without adjusting any of the param- 
eters, see [27] and references therein. In particular 
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FIG. 1: Dynamical quark mass function using the rainbow- 
ladder truncation of Ref. 1221. 



Ftt(Q 2 ) [1, Pm . [Tol [ll[ and the p-n-j, uj-tt-j, and n-jj 
form factors [2J] are well described by this model. We 
therefore expect the model to describe the electromag- 
netic form factors of light vector mesons quite accurately 
as well. 

The corresponding quark propagator functions are 
shown in Fig. [TJ These predictions for the quark mass 
function have been semi-quantitatively confirmed in re- 
cent lattice simulations of QCD [H, [H, |4(3] • Point- wise 
agreement for a range of quark masses requires this in- 
teraction to be flavor-dependent [4l[, and dressing the 



TABLE I: DSE results [22j,|32| for the pseudoscalar and vector 
meson masses and decay constants, together with experimen- 
tal data from Ref. [33l |. unless indicated otherwise. All entries 
are in GeV. 





experiment 


calculated 




(estimates) 


( f fitted) 


m u/d (C = 2 GeV) 


0.003 to 0.006 


0.005 


m s (C = 2 GeV) 


0.095(25) 


0.118 


m c (C, = m c ) 


1.25(9) 


1.30 




0.135, 0.140 


0.138 f 


U 


0.131 


0.13F 


M K 


0.496 


0.497"^ 


fK 


0.160 


0.155 


M p , M w 


0.776, 0.783 


0.742 


fp> 


0.221(2), 0.195(4) 


0.207 


M K * 


0.892 


0.936 


Jk< 


0.224(11) 


0.241 




1.020 


1.074 


u 


0.229(4) 


0.259 




2.980 


2.91 


fnc 


0.335(75) [35] 


0.38 




3.097 




f.J/<S 


0.416(6) 


0.42 
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FIG. 2: Pseudoscalar and vector meson masses and lep- 
tonic decay constants as function of the current quark mass 
(normalized to the physical up/down quark mass). Vertical 
dashed lines indicate the up/down, strange, and charm quark 
masses. 

quark-gluon vertex T l v (g, p) ensures this dependence. The 
consequences of a dressed vertex for the meson BSEs are 
also currently being explored and indications are that 
in the pseudoscalar and vector channels, the effects are 
small [29l . I42I I43I ] . The non-trivial infrared structure of 
the quark-gluon vertex is under investigation using both 
lattice simulations [H, [45|, HU and non-perturbative DSE 
methods [13, EH, but these studies are not yet conclusive. 

B. Quark-mass dependence 



ization point £ = 19 GeV. For the mass region we are 
interested in here, from the chiral limit up to the charm 
quark mass, the fit parameters are 

M PS : C = 0, d = 5.49, C 2 = 5.77, (12) 
My : Co = 0.53, Ci = 6.93, C 2 = 4.88 , (13) 

as illustrated in Fig. [5] For larger quark masses, the fit 
parameters for the pseudoscalar and vector mesons be- 
come closer and closer to each other: a fit on the domain 
m c < m q < 2 nib gives Co = —1.3 and C\ = 7.92 for the 
pseudoscalar mesons, and Co = —0.2 and C\ — 8.02 for 
the vector mesons, with a common parameter C2 = 4.46. 
This reflects the fact that in the heavy-quark limit, the 
pseudoscalar and vector mesons become degenerate: in 
the limit m q — >■ 00 the above fit gives My — Mps — > 
\{CY -Cf)Js/U- 2 *Q. 

The leptonic decay constants increase with the current 
quark mass, both for the pseudoscalar and for the vector 
mesons. Based on the experimental partial decay width 
of vector mesons, it was conjectured [5(| that the vec- 
tor meson decay constants increase with quark mass as 
fy oc \/My. On the other hand [5l|, Coulomb-potential 
models typically give fy oc My, whereas a linear (con- 
fining) potential produces fy ~ constant. Our numerical 
results suggest that both fy and fpg increase approxi- 
mately linearly with quark mass, at least for masses in the 
m c to twice m& range (52j. This is Coulomb-potential- 
like behavior, which may be natural since the effective 
interaction, Eq. [TU1 reduces to one-gluon exchange in the 
ultraviolet region. However, further investigations are 
needed in order to determine the true asymptotic behav- 
ior of the decay constants. 



In recent years, this model has been extended [32j, |49j 
to the charm and bottom quarks as well, in an attempt 
to describe both the light mesons and heavy quarkonia 
within one framework. In Fig. [2] we show our results for 
the vector meson mass and decay constant as function of 
the current quark mass; for comparison we also include 
the evolution of the corresponding pseudoscalar meson 
mass and decay constant. In this and subsequent figures, 
we normalize the current quark masses by the physical 
up and down quark masses of our model; for the strange 
quark we have m s 23 ra u u. 

The charm quark mass m c — 0.827 GeV at the renor- 
malization point £ = 19 GeV is fixed by the experimental 
value of J/ty mass. Again, using one-loop evolution this 
corresponds to m c — 1.30 GeV at £ = m c , see Table Q] 
The decay constant for the J/^S is in good agreement with 
data, and the resulting post-dictions for the rj c mass and 
decay constant are in reasonable agreement with the data 
as well. 

On a limited domain, both the vector and the pseu- 
doscalar meson masses can be fitted reasonably well by 

A^mcson = C + Ci m, + C 2 Ul 2 q , (11) 

where m q is the current quark mass at our renormal- 



C. Frame independence 

The BSE is usually solved in the rest-frame of the me- 
son. However, the calculation of electromagnetic form 
factors in any reference frame entails either the initial 
meson, or the final meson, or both, to have non-zero 3- 
momentum. In a method that is not Poincare covariant 
the wave functions for the moving meson would have to 
be boosted. One of the advantages of the DSE approach 
to hadron physics is its manifest Poincare covariance. 

As an explicit demonstration, we calculate the static 
7T and p properties in a moving frame 2 = (q, 0, 0, i E) 
where q is the 3-momentum of the moving meson [49| . 
Within this frame we solve again the homogeneous BSE, 
Eq. (H|) , and calculate the corresponding electroweak de- 
cay constant. Numerically this is a demanding task, since 
the Lorentz scalar functions of Eq. Q are now functions 
of a radial variable k 2 and two angles 

k-P = i k Ecos a + k q sin a cos (3 , (14) 



2 In the Euclidean metric that we are using here, the rest-frame is 
characterized by = (0, 0,0, iM). 
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FIG. 3: Pion and p mass and decay constant calculated in a 
moving frame, as function of the meson 3-momentum (Figure 
adapted from Ref. [1^]). 



and the integral equation has to be solved in the three in- 
dependent variables k 2 , a, and j3. With current computer 
resources, this can be done without further approxima- 
tions, and the results, shown in Fig. [H are indeed inde- 
pendent of the meson 3-momentum, illustrating that this 
approach is indeed Poincare covariant. We can now use 
this same approach to calculate meson form factors in an 
explicitly covariant manner. 



IV. ELECTROMAGNETIC FORM FACTORS 

A. Vector meson form factors 

Consider the electromagnetic current of a vector meson 
with incoming meson momentum P~ = P p — \Q pi out- 
going meson momentum = P a + \Qa, and incoming 
photon momentum Q M . With this notation, the general 
form for the coup ling of a photon to a vector meson can 
be written as [13, [5311 



A M , PCT (P,Q) 
T^{P,Q) 



3=1 

2p p vTjp-)vUp+) , 



Q 2 

2 M 2 



(15) 
(16) 



v T pAp + ) 



PtS^ \ VUP 1.C17I 



Pu 



p_ 

M 2 



Qp-p, 

X ( Qa + P' a 



" 2 M 2 
Q 2 



p 2 M 2 

Q 2 
2 M 2 



where 



PlM 



(18) 



(19) 



is the transverse projector. The vector meson is on-shell: 
(P~) 2 = (P + ) 2 = -M 2 , where M is the mass of the 
vector meson, and thus P 2 + jQ 2 — —M 2 and P ■ Q = 0. 
This coupling obeys the following relations 



P+A^ a (p, P ') = 0, 
P-A wa {p,p') = 0, 

Qp.Ap,pcr(P,P) = . 



(20) 

(21) 
(22) 



The first two equations simply reflect that the (massive) 
vector mesons are transverse; the last equation follows 
from current conservation. 

The electric, magnetic, and quadrupole form factors 
Ge , Gm , Gq can be expressed in terms of these scalar 
functions F; 



2 Q 2 



3 4M 2 



G E (Q 2 ) = Fi(Q 2 )- 
G M (Q 2 ) = -F2(Q 2 ) , 
Gq(Q 2 ) = F,(Q 2 ) + F 2 (Q 2 ) 



Gq(Q 2 



(23) 
(24) 



Q 2 



4M 2 



F 3 (Q 2 ) 



(25) 



The electric monopole moment (i.e. the electric charge), 
magnetic dipole moment and the electric quadrupole mo- 
ment follow from the values of these form factors in the 
limit Q 2 -> 

(26) 
(27) 
(28) 

Here the magnetic moment \x and the quadrupole mo- 
ment Q are introduced; the electric charge is 1 in terms 
of the fundamental charge unit e. For point-like vec- 
tor particles, the magnetic and quadrupole moments are 



G E (Q 2 


= 0) 


= 1 , 


Gm(Q 2 


= 0) 


= V , 


Gq(Q 2 


= 0) 


= Q 



[i = 2 in units of e/2My and Q 
respectively (54| . 



T in units of e/M, 



v 



B. Impulse approximation 

The generalized impulse approximation allows electro- 
magnetic processes to be described in terms of dressed 
quark propagators, bound state BSAs, and the dressed 
17177-vertex, see Fig. 01 In combination with ladder- 
rainbow truncation for the vertices and the quark prop- 
agators, it ensures electromagnetic current conserva- 
tion [1,0], see also Sec. lIVCl below. Phenomenologically, 
this approximation has proved to be very successful in 
describing the pion electromagnetic form factor [l(| [H| • 

Consider for example the 3-point function describing 
the coupling of a photon with momentum Q to a vector 
meson ab, with initial and final momenta P ± Q/2. This 
interaction can be written as the sum of two terms 



A 



ab 



r (P,Q) 



Qa a aab 



Q b K b 



(29) 
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FIG. 4: Meson form factor in impulse approximation 



where Q is the quark or antiquark electric charge, and 
where A aha (P, Q) and A ahh (P, Q) describe the coupling of 
a photon to the quark (a) and antiquark (b) respectively. 
In impulse approximation, these couplings are given by 

A^V(P,Q) = iN c [Tr[ri(q^q + ) X f(q + ,q) 

Jk 

x^})" 1 ^,?-)], (30) 

with q = k - P/2 and q± = k + P/2 ± Q/2, and simi- 
larly for Af> pa . 

For the p mesons it is sufficient to calculate the cou- 
pling of the photon to a single quark, a direct conse- 
quence of isospin invariance, while for the K* mesons we 
add contributions from photon coupling to the quark and 
the antiquark. Thus for the charged K* the form factors 
are given by 

F.r + (Q 2 ) = lFr s {Q 2 )~\Fr s {Q 2 ) (31) 
= |f™«(Q2) + \F? sn {Q 2 ) , (32) 

and similarly for the neutral K* we have 

Fr"(Q 2 ) = -\Fr s {Q 2 ) + \Fr{Q 2 ). (33) 



C. Quark-photon vertex 



truncation [f|0] 

r M (pout,Pin;Q) = ^2 7^+ / K(p out ,p in ; 

^out j Mn) 

Jk 

Q) S(k in ) , (35) 

with p ou t and pi n the outgoing and incoming quark mo- 
menta, respectively, and similarly for fc out and k mi with 
Pout — Pin = fcout — fcin = Q- The kernel K is the same 
kernel as used in the meson BSE, defined under Eq. 

Note that solutions of the homogeneous version of 
Eq. (|35p define qq vector meson bound states with masses 
My = —Q 2 at discrete timelike momenta Q 2 . It follows 
that r M has poles at those locations and, in the neigh- 
borhood of Q 2 = —M v , behaves like @ 

P , x rV( Pout ,p in ) /ytfy 

r M bout, Pin) q 2 + m 2 ' *■ ' 

where TY is the qq vector meson BSA, and /v the cor- 
responding electroweak decay constant. The fact that 
the dressed qqj- vertex exhibits these vector meson poles 
explains the success of naive vector-meson-dominance 
[VMD] models; the effects of intermediate vector meson 
states on electromagnetic processes can be unambigu- 
ously incorporated by using the properly dressed qqj- 
vertex rather than the bare vertex 7^ [|| . 

V. DISCUSSION OF NUMERICAL RESULTS 

Numerical solutions for the quark propagator, vec- 
tor meson BSAs and dressed quark-gluon vertex can 
now be used in Eq. (|30f to calculate the electromag- 
netic form factors. We explicitly solve the respective 
(in) homogeneous BSEs for the meson BSAs and for the 
gg7-vertex in the corresponding momentum frame, thus 
avoiding any need for interpolation or extrapolation of 
the numerical solutions of the BSEs. This does mean 
that we have to solve the meson BSE for each value of 
Q 2 , but the advantages of not having to extrapolate the 
numerical BSAs outweighs the additional numerical ef- 
fort of repeatedly solving the meson BSE. Further details 
are given in the appendix. 



The impulse approximation and rainbow-ladder trun- 
cation together lead to current conservation only if 
the quark-photon vertex T M satisfies the vector Ward- 
Takahashi identity 0] 



iQ M r„(fc + ,fc_;Q) = S-^k+^-S-^k 



\Q). 
(34) 



This identity can be satisfied by the Ball-Chiu 
Ansatz (5f| for the quark-photon vertex; a more consis- 
tent approach, which we follow, is to use the solution of 
the inhomogeneous BSE for the gyy-vertex in the ladder 



A. Rho and K* form factors 

Results for the charge radius, 

dG E (Q 2 ) 



= -6 



dQ 2 



(37) 



Q 2 =0 



the magnetic moment /1 = Gm(0), and the magnetic 
quadrupole moment Q — Gq(0), are presented in Ta- 
ble [TTJ Experimental data on these static properties is 
absent, and so results from other theoretical models are 
included in the table. 
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TABLE II: Results for the p and K* meson charge radii (r v ) 
in fm 2 , magnetic moments, and quadrupole moments, com- 
pared to other calculations. For comparison, we also include 
the results for the pseudoscalar charge radii, calculated within 
the same model, where available, as well as the experimental 
pseudoscalar charge radii. 



7r, p meson 


r 2 


r 2 




Q 


current DSE [6, 32] 


0.44 


0.54 


2.01 


-0.41 


previous DSE [4, 15} 


0.31 


0.37 


2.69 


-0.84 


Covariant QM [13] 




0.37 


2.14 


-0.79 


Lightcone QM [5, 17} 


0.43 


0.27 


1.92 


-0.43 


experimental [33] 


0.452(11) 








us meson 


r 2 K 


r 2 K * 


A* 


Q 


current DSE calc. 


0.38 


0.43 


2.23 


-0.38 


previous DSE [15] 


0.28 


0.29 


2.37 


-0.62 


experimental [33] 


0.314(35) 








ds meson 


r 2 


r 2 




Q 


current DSE calc. 


-0.09 


-0.08 


-0.26 


0.01 


previous DSE [15] 


-0.03 


-0.05 


-0.40 


0.11 


experimental [33] 


-0.077(10) 








cc meson 


r 2 


2 

r J/i> 




Q 


current DSE [32] 


0.048(4) 


0.052(3) 


2.13(4) 


-0.28(1) 


lattice [20] 


0.063(1) 


0.066(2) 


2.10(3) 


-0.23(2) 



Our calculations give a significantly larger value for 
the charge radius of the p and the charged K* mesons 
compared to a previous DSE calculation lfj, see Ta- 
ble HU This is the case not only for the vector mesons, 
but also for the pion and kaon charge radii. This dif- 
ference can largely be attributed @ to the fact that we 
employ a dressed quark-photon vertex that has poles in 
the neighborhood of Q 2 — -My. Thus the effects from 
intermediate vector meson states are unambiguously in- 
cluded in our calculation. On the other hand, Ref (l5| 
uses the Ball-Chiu Ansatz for the gg7-vertex. Since our 
results compare favorably with the experimental charge 
radii, we expect our results for the vector radii to be more 
realistic than those of Ref. [HI . 

The quark model [QM] calculations of Refs. [l3| and 
[ijj report a considerably smaller value for (r 2 ). Again, 
this is in part due to the fact that neither of these cal- 
culations incorporate the effects of vector meson poles in 
the quark-photon vertex [||. 

Another significant difference between our results and 
those of Ref. [15[ is that the latter predicts a zero cross- 
ing of Ge(Q 2 ) at about Q 2 « 1.7 GeV 2 for the p-meson, 
whereas we do not find any evidence for such a zero 
crossing below Q 2 = 2.5 GeV 2 . Extrapolating our nu- 
merical results for Ge(Q 2 ) using a (2,3)-Pade fit suggests 
a zero crossing around Q 2 « 3.8 GeV 2 . For comparison, 
Ref. 13 1 predicts a zero crossing at Q 2 » 2.9 GeV 2 , and 
Ref. 17 1 at Q 2 ~ 5.5 GeV 2 , whereas neither of the sum 
rule analyses [l8| and Ref. EH predict a zero crossing of 
Ge(Q 2 ) below Q 2 = 5 GeV . Based on our calculations, 
and on the rather wide range of predictions from other 



TABLE III: Our results for the p meson electric, magnetic, 
and quadrupole form factors Ge,m,q(Q 2 ) at Q 2 = 1 GeV 2 
and at Q 2 = 2 GeV 2 , compared to previous DSE [15|], light- 
cone 171, and sum rule calculations Hal. 





Q 2 


= 1 GeV 


2 


Q 2 


= 2 GeV 


_> 




Ge 


Gm 


Gq 


Ge 


Gm 


Gq 


current 


0.22 


0.57 


-0.11 


0.08 


0.27 


-0.05 


Ref. [15] 


0.17 


0.85 


-0.51 


-0.04 


0.45 


-0.32 


Ref. [17] 


0.38 


0.93 


0.23 


0.18 


0.59 


-0.15 


Ref. [19] 


0.25 


0.58 


0.49 


0.13 


0.28 


-0.24 



calculations, we conclude that it is unlikely for Ge(Q 2 ) 
to have a zero crossing below Q 2 ss 3 GeV 2 . 

The values we obtain for the magnetic and quadrupole 
moments of the p-meson, p and Q, are very similar to 
those of Ref. |17| . and significantly smaller than those of 
[IH and Refs. |15|. However, the Q 2 evolution of both 
Gm and Gq is quite different than that of Ref. [ItJ (see 
Table llllj) , again most likely due to the fact that VMD 
effects are not properly accounted for in Ref. [l7j |. 

A recent sum rule analysis [l6[ obtained p = 2.0 ± 0.3 
for the magnetic moment, which, given the large error 
bars, is not accurate enough to discriminate between dif- 
ferent calculations. The sum rule analysis of p form fac- 
tors at Q 2 = 1 GeV 2 and 2 GeV 2 ^ seems to support 
our calculation, see Table IrTrl at least for Ge and Gm- 
Our results for the quadrupole form factor however are 
almost a factor of five smaller than those of Ref. [lj|. 
Clearly the quadrupole form factor is most sensitive to 
the details of the dynamics. 

As we can see in Figs. [5] and HI not only Ge(Q 2 ), but 
also Gm(<2 2 ) and Gq(Q 2 ) have poles in the timelike re- 
gion as one approaches Q 2 = —M 2 ; in addition, the K* 
form factors also have poles at Q 2 = — Mj. In Fig. [S] 
we also show the electromagnetic form factors of a ficti- 
tious ss-like vector meson with the photon coupled (with 
charge one) to only the s-quark, but not to the s-quark 3 . 
The form factors of this fictitious ss-like vector meson 
have poles at Q 2 = —M^. For comparison we also show 
a pure VMD model for the electric form factor 



G E (Q 2 



M 2 



M 2 



Q 2 



(38) 



for both the p and the ss-like meson in the top panel 
of Fig. [5] Such a VMD model works remarkably well 
for the pion electromagnetic form factor, at least up to 
spacelike Q 2 values of about 4 GeV 2 . Fig. [5] shows that 
our form factors of the vector mesons deviate significantly 
from a simple VMD curve. Only in the timelike region, 
near the actual vector meson pole, does the VMD curve 



3 The form factors of the physical <j> meson are trivially zero, as 
are the form factors of the neutral p meson. 
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Q 2 [GeV 2 ] Q 2 [GeV 2 ] 
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Q 2 [GeV 2 ] 
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-0.2 0.2 0.4 0.6 0.8 1 
Q 2 [GeV 2 ] 




Q 2 [GeV 2 ] Q 2 [GeV 2 ] 



FIG. 5: Numerical results for G E (Q 2 ) (top), G M (Q 2 ) (mid- 
dle), and Gq(Q 2 ) (bottom) for the p meson and a (fictitious) 
ss state. For comparison, we also show the VMD result 
(dashed and dot-dashed curves) for Ge(Q 2 )- 



resemble our results. In the spacelike region, our results 
drop significantly faster than a VMD form factor. 

Our results for the charged K* meson form factors are 
qualitatively similar to those for the p meson. The neu- 
tral K* form factors are most sensitive to details of the 
calculation, since it depends on the difference between 



FIG. 6: Numerical results for the neutral and charged K* 
form factors Ge(Q 2 ) (top), Gm(Q 2 ) (middle), and Gq(Q 2 ) 
(bottom). 

the up/down quarks and the strange quarks. Therefore 
we expect these form factors to show more model depen- 
dence than the presented form factors of the charged K* 
and p mesons, and this is indeed what we see in Table [TT1 
again in particular for the quadrupolc moment. 

Finally, from Table |TT] it is interesting to note that all 
studies, with the exception of the light-cone QM calcu- 
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lation of Ref . [17| , find the charge radii of (charged) vec- 
tor mesons to be larger than those of the correspond- 
ing pseudoscalar mesons. This trend was recently con- 
firmed in lattice calculations for light quarks [21| and 
for charmonium-like states [20J, with the photon cou- 
pled to the quark only, not to the antiquark, of the 
charmonium state. Also a recent non-relativistic QM 
calculation of such charmonium-like states gives a vec- 
tor charge radius that is larger than the pseudoscalar 
charge radius [56[ . This means that the vector states are 
broader than the corresponding pseudoscalar states, as- 
suming that the charge distribution is indicative of the 
physical size of the bound state. This agrees with the 
naive intuition that a more tightly bound state is more 
compact than a heavier state with the same constituents. 



"3 1.5 

o 



o 
O 



0.5 



i 



x electric charge 

♦ magnetic moment 

• quadrupole moment 



;:«o<xxxx x x x 



50 100 
m / m , . 
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150 200 
[dimensionless] 



250 



B. Quark-mass dependence 

In addition to the electromagnetic form factors of the 
physical p meson, we have also calculated the quark- 
mass dependence of these form factors. For simplicity 
we restrict ourselves to equal-mass mesons, i.e. qq bound 
states, and for the electromagnetic properties we couple 
the photon to only the quark q (with charge one) , not to 
the antiquark q. We refer to these form factors as the 
single-quark form factors; even though they are unphysi- 
cal, they are well-defined and allow for comparisons with 
other theoretical and computational studies of the vector 
meson form factors. 

The static electromagnetic properties are plotted in 
Fig. [JJ as function of the current quark mass, normal- 
ized by the up/down current quark mass; for the strange 
quark we have m s « 23 m u u and for the charm mass 
m c « 224 m u / d . Our numerical errors grow with in- 
creasing meson mass (i.e. increasing quark mass), mainly 
because the momentum p 2 of the quark propagator in 
the Bethe-Salpeter integrals and in the triangle diagram 
spans an increasingly large domain in the complex p 2 - 
plane with increasing meson mass. The analytic contin- 
uation of the solution of the quark DSE, Eq. fT]), from 
the spacelike (Euclidean) axis to this complex momen- 
tum domain becomes numerically cumbersome and inac- 
curate for large meson masses. In our calculations this 
is reflected by the fact that for the charm quark we have 
an estimated 2% to 4% numerical uncertainty in the mo- 
ments, and a numerical uncertainty of 6% to 8% in the 
radii. 

The magnetic moment (i turns out to be almost in- 
dependent of the current quark mass, Fig. [JJ For < 
m q < 4 m s rs 100 m u / d the magnetic moment increases 
with quark mass from \x = 2.01 to p = 2.12; above 4 m s 
the numerical uncertainty starts to increase and the mag- 
netic moment is basically independent of the quark mass 
within our numerical error bars of a few percent. 

This mass dependence can naively be understood by 
considering a simple constituent quark model, in which a 
vector meson is a bound state of two quarks in an S-wave 




x x charge radius 
♦ ♦ magnetic radius 
- - VMD curve 



100 
m /m ,, 

q up/down 



150 200 
[dimensionless] 



250 



FIG. 7: Numerical results for the magnetic and quadrupole 
moments (dimensionless, top) and for the charge radius (in 
fermi, bottom) as function of the current quark mass. Vertical 
dashed lines indicate the up/down, strange, and charm quark 
masses. 



with the spins aligned. In such a model one expects the 
magnetic moment to be proportional to that of its con- 
stituents, i.e. nv oc fig, where (j, q is the magnetic moment 
of the constituent quark. As long as the quark magnetic 
moment is only weakly dependent on the current quark 
mass, so is the magnetic moment of the bound state. 

Since both the charge, Ge(0), and the magnetic mo- 
ment, Gm(0), are (almost) independent of the quark 
mass, one might expect that also the quadrupole mo- 
ment, Gq(0), depends only weakly on the quark mass. 
However, our calculation shows that this is not the case: 
The quadrupole moment Q decreases monotonically with 
rn 



and is reduced by about 25% at m q 



50 m u / d . Above this quark mass, the quadrupole mo- 
ment continues to decrease with quark mass, but at a 
slower rate. The mass dependence of the quadrupole 
moment may not be that surprising on realizing that our 
results, —0.41 < Q < —0.27 depending on quark mass, 
deviate significantly from the canonical value, Q = — 1. 
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Furthermore, of all the static electromagnetic properties, 
the quadrupole moment shows the largest model depen- 
dence, see Table QJ indicating that the details of the 
dynamics are important for this quantity. 

The fact that we find a nontrivial quadrupole moment 
for all quark masses, indicates that the naive picture of 
a vector meson as a nonrelativistic bound state of two 
quarks in an S-wave is too simple. Clearly, there is a sig- 
nificant amount of quark orbital momentum in the vector 
mesons, which in our approach is incorporated in the me- 
son BSAs. This is related to the fact that we use method 
that is explicitly frame- independent: orbital angular mo- 
mentum is not a Poincare invariant quantity. 

In the bottom panel of Fig. [7] we show the mass de- 
pendence of the charge and magnetic radii (defined anal- 
ogously) of the vector mesons. Just like the pion charge 
radius [49|], both the charge and the magnetic radius de- 
creases with increasing quark mass. Our results for (r 2 ) 3 
are qualitatively similar to a VMD curve, and over the 
entire mass range the charge radius can in fact be rea- 
sonably well described by a VMD curve with a constant 
shift of about 0.07 fm. Of course that means that the 
relative deviation from a VMD model increases with in- 
creasing quark mass, again just as in the case of the pion 
form factor. 

Finally, we consider the single-quark transition form 
factors of the rj c and J/*f? mesons, at m c w 224?7j n / C 2- 
This allows us to compare our results with the lattice 
simulations of Ref. [2(| • Both the magnetic moment and 
the quadrupole moment are in good agreement with the 
lattice data, see the bottom row of Table [H] Also the 
charge radii are in reasonable agreement, with the vector 
state being slightly larger than the pseudoscalar state. 

In this paper we are interested in the quark-core con- 
tribution to mesonic observables, and do not include pion 
loop effects. However, the form factors, and their quark- 
mass dependence, are sensitive to pion loops, in partic- 
ular at small quark masses. Pion loops lead to correc- 
tions to the physical p and ir charge radii of about 10% 
to 15% [57, 58], but their contributions decrease rapidly 
with increasing quark mass. Similarly, we expect that 
at the physical up/down quark masses, the other form 
factors, Gm and Gq will also receive corrections of the 
order of 10% due to pion loops. That means that the 
mass dependence of the magnetic moment will be domi- 
nated by the pion loop corrections, given the very weak 
quark-mass dependence we find here. However, we ex- 
pect these corrections to become negligible for masses 
around m s and above. 



VI. CONCLUSIONS 

We have calculated the electromagnetic form factors 
of the p meson, and of both the charged and the neutral 
K* mesons. Our method is explicitly Poincare invariant, 
and we have demonstrated explicitly that physical ob- 
servables are frame- independent. By dressing the quark- 



photon vertex we guarantee electromagnetic current con- 
servations; furthermore, our dressed quark-photon vertex 
exhibits poles in the timelike region, corresponding to in- 
termediate vector mesons. Exactly the same method, us- 
ing the same model for the effective interaction [22j, has 
been used quite successfully to describe the pion elec- 
tromagnetic form factor 0, [HJ and a plethora of other 
light meson observables within about 10% to 15% [27| . 
We therefore expect our results to have a similar level of 
accuracy. 

Compared to other calculations, we find a stronger Q 2 - 
dependence, mostly due to the fact that we have incor- 
porated unambiguously VMD effects in the quark-photon 
vertex, which have been neglected in Refs. [H, [H, E3 • 
Our results favor a magnetic moment close to two, 
in agreement with sum rule analysis (T5 , Il9j . and a 
quadrupole moment Q ~ —0.4, which is similar to the re- 
sult of Ref. [l7| , but significantly smaller than suggested 
by other calculations [lj, [jja, [lj|. Our results for the 
quadrupole moment indicate that there is a significant 
amount of quark orbital angular momentum in the vec- 
tor mesons. 

The magnetic moment is almost independent of the 
quark mass for mesons of equal-mass constituents. The 
quadrupole moment decreases with increasing quark 
mass. Also the shape of the form factors changes: since 
the VMD pole in the timclikc region shifts further away 
from Q 2 = with increasing quark mass, the form factors 
become less steep, and the radii decrease with increasing 
quark mass. Over the entire quark mass range, from the 
chiral limit to the charm quark mass, the charge radius 
can be reasonably described by a VMD curve with a con- 
stant shift of about 0.07 fm. A similar behavior of the 
form factors, and of their quark-mass dep endence, has 
been obtained using the model of Ref. [59( for the effec- 
tive quark-antiquark interaction. 

Ideally experiments would guide us in order to discrim- 
inate between different models and different calculation 
methods, but it is unlikely that these form factors can 
be measured in the near future. Lacking reliable experi- 
mental input, it would be very useful to have (quenched) 
lattice data for the vector meson form factors at light 
quark masses. Currently, the only accurate lattice data 
available are at the charm mass [2fJ ; these lattice data are 
in reasonable agreement with our results at m c . Accu- 
rate lattice simulations at light quark masses are needed 
in order to discriminate between different model calcula- 
tions. A beginning has been made in Ref. [2l|, but the 
error bars have to be reduced significantly in order to 
make a detailed comparison meaningful. 
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APPENDIX A: NUMERICAL APPROACH 

The form factors are calculated in impulse approxima- 
tion 

xS^kp)- 1 ^,*-)], (Al) 

with k P = k-P/2 and k± = k + P/2 ± Q/2. That 
means that we need numerical solutions of the rainbow- 
ladder BSE for 

• quark-photon vertex 

• incoming vector meson 

• outgoing vector meson 

in addition to the solution of the quark DSE in rainbow 
truncation. If we choose the momentum routing and in- 
tegration variables carefully, we can arrange the integra- 
tion grids such that we do not need any interpolation or 
extrapolation in the the final (triangle) loop integral for 
the form factor. That significantly reduces numerical er- 
rors, in particular possible systematic errors introduced 
by the extrapolation. 

We use the momentum frame 

P^ = (0,0,0,P), (A2) 
Qti = (Q, 0,0,0) , (A3) 

with P 2 + jQ 2 — —M 2 . Depending on the value of Q 2 , 
either P or Q or both are imaginary. In addition, we 
have the integration momentum 

= /c(cos(0),sin(0)cos(0),sin(0)sin(0),O) , (A4) 

and corresponding measure 

r d 4 k f°° k 3 dk r . 3 m r . ^ 

(A5) 

We use Eqs. (fl"6l) -(fl"8 |) for the general structure of the 
coupling of a photon to an on-shell vector meson, and 
perform the traces analytically to obtain expressions for 
the form factors Fi{Q 2 ). Subsequently, we use Eqs. (|23|) - 
(|2l))) to convert the functions Ft to the more conventional 
electric, magnetic, and quadrupole form factors Ge, C?m, 
and Gq. 



1. Quark-photon vertex 

The gg7-vertex is the solution of an inhomogeneous 
BSE, which in ladder truncation we can write as 

r M (fc_,fe + ;Q) = Z 2 ^ + ^J-0^4na{(k- q ) 2 ) 

x D pa {k - q) 7 P Xn(q-,q+;Q) la , 

(A6) 

where 

Xti(q-,q+;Q) = S(q-) r M (g_,g + ;Q) S(q+) , (A7) 

with k± = k + P/2 ± Q/2, and similarly for q±. Both k 
and q are real Euclidean vectors, and the incoming and 
outgoing quarks are always of the same flavor, so we can 
drop any flavor indices. 

In order to solve such a BSE, we decompose the func- 
tion x into its Dirac components, and project out a set of 
coupled integral equations for its scalar component func- 
tions Fj. (Equivalently, one could decompose T into its 
Dirac components, and solve the coupled integral equa- 
tions for its components, but it turns out that solving 
the BSE for \ is 1S roughly factor of three faster than 
solving the BSE for T.) We solve these coupled integral 
equations by iteration, after discretizing the angular and 
radial variables. 

The most general form of the quark-photon vertex re- 
quires twelve Dirac structures. Four of these covariants 
represent the longitudinal components which are com- 
pletely specified by the Ward-Takahashi identity in terms 
of the (inverse) quark propagator and they do not con- 
tribute to elastic form factors. The transverse part of the 
vertex x can be decomposed into eight components 

X^(k-,k + ;Q) = 

8 

Y,Ti(k + ±P,Q)F i {k 2 ,6,<f>) , (A8) 
»=i 

where the functions F{ depend on two angles, because we 
solve it in exactly the same frame and variables as it is 
used in the loop integral for the form factor; alternatively, 
we could write it as a function of fc 2 , k ■ P, and k ■ Q. 

2. Vector meson BSAs 

The incoming vector meson, with momentum P — 
and flavor labels ab, is the solution of the homogeneous 
BSE 

S-\k+) xf(k+,k P ;P-lQ) S^(k P ) - 

xj aX f(q+,qp;P-kQ)^, (A9) 
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with k+ = k + P/2 + Q/2, k P = k - P/2, and similarly 
for q + and qp. The vector meson is on its mass-shell: 

p2 + 1 Q 2 = _ M 2_ 

Again, we decompose the transverse vertex function 
X into eight components, but use a slightly different de- 
composition than for the quark-photon vertex, because 
the momentum arguments are different. For the vector 
meson BSAs we use 



Xp {k + ,k P ;P-\Q) = 

8 

5^T j (fc P ,P-!Q).Fi(fc 2 ,0,, 



(A10) 



i.e. we use the momentum partitioning where the to- 
tal incoming meson momentum flows into the outgoing 
quark leg; or in other words, we use the incoming quark 
momentum as the "relative momentum" in the decompo- 
sition (t) = 1), though not as the integration momentum. 
And again, the functions Fi depend on two angles. As 
a check on our numerics, we calculate the leptonic decay 
constant in this frame as well, which gives us an indica- 
tion of the numerical errors in the BSAs. 



We do not explicitly use any algebraic relation between 
X and x 



x(p,p) = Ic- 1 x (-p,-p)c 



transpose 



(All) 



because the arguments of x an d X are different. Instead, 
we simply solve both for x and for x m the appropriate 
frame. That is, we calculate the BSA of the outgoing 
vector meson as the solution of 



Sb\k P ) xl a {k Pl k_-P+\Q)S-\k^) 
4 f d 4 q 



Ana({k~q) 2 ) D aP {k-q) 



(2tt)< 

x 7a x b p{qp,q-;P + \Q) m 



(A12) 



again with fc_ = k + P/2 — Q/2, kp = k — P/2, and sim- 
ilarly for g_ and qp. We solve this equation for x m 
basically the same manner as the BSE for x- This way 
we avoid the need for interpolation and extrapolation on 
the vertex functions in the triangle diagram for the elec- 
tromagnetic form factors. 
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